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Abstract 

We consider pyramids made of one-dimensional pieces of fixed integer length a and which may 
have pairwise overlaps of integer length from 1 to a. We prove that the number of pyramids of 
size m, i.e. consisting of m pieces, equals ( a ™Ii) f° r ea ch a > 2. This generalises a well known 
result for a = 2. A bijective correspondence between so-called right (or left) pyramids and 
a-ary trees is pointed out, and it is shown that asymptotically the average width of pyramids is 
proportional to the square root of the size. 
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1 Introduction 



Solutions of the enumeration problem for a variety of lattice animals have been obtained in 
recent years, although the problem concerning general animals on the square or the triangular 
lattice both remain unsolved. The problems considered in this area are interesting in their own 
right from a combinatorial point of view, some of them being equivalent to other well known 
combinatorial problems, and frequently they are inspired by concrete problems in other fields. 
Thus, for the standard model of site percolation on a lattice [13], the connected percolation 
clusters are lattice animals and their combinatorial properties of imminent importance for the 
critical properties of the model [9] . Viewing the points in a lattice as centres of the elementary 
cells of the dual lattice, an animal can be identified with an edge-connected set of elementary cells 
on the dual lattice, also called a polyomino. Thus, enumerating lattice animals and polyominoes 
amounts to the same problem, although the motivation for studying a certain class of animals 
and its polyomino counterpart may be quite different. Early enumeration results for polyominoes 
can be found in [16,23]. For more recent results see [6,26] and references given there. 

A particular, much studied, problem is that of directed animals on a square lattice, first solved 
by Dhar in [10], and later by a number of authors using different methods [3-6, 14, 15, 19, 20]. 
A directed animal on the square lattice is a set of points on the lattice such that any point in 
the set is the end point of a lattice path starting at the origin all of whose steps are directed 
towards either east or north and all of whose points are contained in the set. The perhaps most 
elegant solution to this problem [6] is obtained by observing [24] that if the lattice is rotated 
through an angle it/ 2 counterclockwise and each point in a directed animal is replaced by a 
suitable dimer, one obtains a pyramid of dimers, the detailed definition of which is given below, 
restricted such that no dimer is placed directly on top of another. The generating function for 
such pyramids can then be obtained rather simply as a solution of an algebraic equation. In 
turn, the solution to the directed animal problem on the triangular lattice can be obtained using 
that those animals correspond to general pyramids of dimers [24] . 

In this paper we address the problem of enumerating pyramids whose pieces are of fixed but 
arbitrary integer length a instead of dimers (which correspond to a = 2). These may, of course, 
be viewed as a particular type of polyominoes on the square lattice, or one may think of them as 
connected, planar LEGOs made of 1 x a-pieces and which are obtained by dropping successively 
pieces from above so that the resulting configuration is connected. 

More precisely, we shall consider heaps in the sense of Viennot [25] whose basic pieces are 
one-dimensional and have fixed integer length a and whose concurrency relation is defined by 
assuming each piece to be an interval ]s, s + a[, s £ Z, and two intervals a, (3 are concurrent if 
and only if a f~l (3 ^ 0. Thus a heap, in this article, can be thought of as being obtained by 
dropping a finite number of pieces towards a horizontal axis. Recall that a heap is a pyramid if 
it has a unique bottom piece. We call a pyramid p a right s -pyramid, if the bottom piece covers 
the interval ]s, s + a[ and is a leftmost piece in p. Similarly, p is a left s-pyramid if the bottom 
piece covers the interval ]s — a, s[ and is a rightmost piece in p. 

When using the term pyramid it will henceforth be assumed, unless otherwise stated, that 
its bottom piece covers the interval ]0, a[. The number of pieces in a pyramid p will be called 
its size and is denoted by \p\ . 

The main result of this note is the following. 

Theorem 1.1. Given a > 2, the number of pyramids of size m equals (^-i) ■ 

The result is well known for a = 2 [6] . We reconsider this case in Section 2 for the purpose 
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of illustrating our method of argument which, in particular, involves establishing a bijective 
correspondence between pyramids of size m and closed walks on the integers of length m. In 
Sections 3 and 4 we generalise in two steps this correspondence to the case a > 3 from which 
the main result will follow. It is worth noting that to obtain this result we do not rely on gen- 
erating function techniques whose applicability seems to be restricted to the dimer case. Those 
techniques, on the other hand, are used to determine the asymptotic behaviour of the average 
width of pyramids of large size in Proposition 4.3. We also point out a bijective correspondence 
between right pyramids of size m and a-ary trees with m nodes. In Section 5 we conclude with 
some numerical results and comments concerning the growth rate of the number of general pla- 
nar LEGOs as a function of size and on the dependence of the exponential growth constant on 
the size of the building blocks. 



2 The dimer case 

In this section we assume a = 2. Hence the pieces in this case can be thought of as dimers. 
We first note the following decomposition property. 

Lemma 2.1. There is a bijective correspondence between pyramids of size m > 1 and sequences 
(Pi)P2>- • ■ >Pr) of pyramids such that pi is a right 0-pyramid if i is odd and a left 1-pyramid if i 
is even, and such that \pi \ + ■ ■ ■ + \p r \ = m. 

Proof. If the pyramid p is not a right 0-pyramid there is a lowest piece in p above the interval 
] — 1 , 1 [ and this piece is the bottom piece of a unique proper sub- pyramid p' and we can write 
p = Pi Op', with notation as in [25], where p% is a right 0-pyramid. If p' is not a left 1-pyramid 
it contains a unique lowest piece above the interval ]0,2[ and we have p' = p2 p" , where 
P2 is a left 1-pyramid. Repeating the argument the claim follows. See also [6] for a similar 
decomposition. □ 

Definition 2.2. A finite sequence of O's and l's will be called a string and by a (n,m) -string 
we mean a string of length n with m l's. A (2m, m)-string X1X2 ■ ■ ■ X2 m is called positive if 

s 

t B = ^(2x„ - 1) 

u=l 

is non-negative for all s = 1, . . . , 2m, i.e. the number of O's in x\ . . . x s at most equals the number 
of l's in x\ . . . x s for each s. 

Note that a positive (2m, m)-string necessarily begins with a 1 and ends with a 0. There is a 
natural correspondence between strings and nearest neighbouring walks on the integers starting 
at where each corresponds to a left-step and each 1 to a right- step. Positive strings then 
correspond to walks on the non-negative integers starting at 0. 

Lemma 2.3. There is a bijective correspondence between positive (2m, m) -strings and right 
0-pyramids. 

Proof. Let p be a right 0-pyramid of size m. We construct inductively the corresponding positive 
string x\ . . . X2 m together with a sequence p^ . . . p( 2m ) of 0-pyramids that are sub-pyramids of 
p such that p( 2m ) = p as follows. 
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Let x\ = 1 and p^ be the bottom piece of p. Assume x\ . . . x s and p^ . . . p^ have been 
constructed. If a piece above the interval ]t s ,t s + 1[ can be dropped onto p^ to obtain a sub- 
pyramid of p we let be that pyramid and set x s +i = 1. Otherwise, set = p^ and 



x s +i = 0. Here t s is given as in Definition 2.2 and one readily checks that at any stage t s is less 
than the width of p( s \ i.e. the length of the projection of p^ onto the horizontal axis, and that 
the size of p^ equals the number of l's in x\ . . . x s . Indeed, by construction, any piece that can 
be dropped onto p^ to obtain a sub-pyramid of p is above some interval contained in ]0,t s + 1[, 
and at any stage we choose the rightmost of those pieces to obtain p( s+1 \ It follows that the 
so obtained sequence x\ . . . x<im after 2m steps is a positive (2m, m)-string since otherwise the 
number of l's would be less than m and t2m would hence be negative, which is not possible. 

If p p', the corresponding sequences p^ 1 ' . . . p( 2m ) and p'^ . . . p'( 2m ) will deviate at some 
minimal step s, 1 < s < 2m, and it follows that the corresponding strings also deviate at step 
s. On the other hand, any positive (2m, m)-string x\ . . . X2m can be obtained by the described 
procedure from the right 0-pyramid p obtained by successively dropping a piece above those 
intervals }t s ,t s + 1[ for which x s = 1, with the convention to = 0. This concludes the proof. □ 

Lemma 2.4. Any (2m, m)-string w = x\ . . . X2 m starting with x\ = 1 can be written in a unique 
way by juxtaposition as 

w = W\ . . . w r 

where Wi is a positive 0-string if % is odd, and w~ l is a positive 0-string if i is even. Here w^ 1 
denotes the string obtained from Wi by reversing its order. 

Proof. Using the correspondence between strings and nearest neighbouring walks on the integers, 
we see that the statement amounts to asserting the obvious unique decomposition of a walk 
starting and ending at into an alternating sequence of walks on the non-negative, respectively 
the non-positive, integers. □ 

We are now in a position to derive the following result which, in particular, proves Theo- 
rem [TTT] in case a = 2. 



Proposition 2.5. There is a bijective correspondence between pyramids of size m and (2m, m)- 
strings starting with 1. In particular, the number of pyramids of size m equals (^Zi)- 

Proof. Since, obviously, there is a bijective correspondence between left and right 0-pyramids 
of given size and since reversal of ordering of a string is injective, the claimed correspondence 
follows from the preceding three lemmas. The last statement follows by noting that a (2m, m)- 
string starting with 1 is uniquely determined by the position of the remaining l's among the 
remaining 2m — 1 entries of the string. □ 



3 Decomposition of pyramids and strings for a > 3 

In this and the subsequent section we assume a is fixed and larger than or equal to 3. Accordingly, 
we extend the notion of positive strings as follows. 

Definition 3.1. An (am, m)-string x\X2 ■ ■ ■ x am is called positive if 

s 

t s = ^(ax u - 1) 

u=l 
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Figure 1: The ten 2-pyramids of size 3 and their description as (6, 3)-strings 

is non-negative for all s = 1, . . . , am. Moreover, a string is called negative if the reversed string 
is positive. 

Note that a positive string necessarily begins with 1 and ends with at least a — 1 consecutive 

O's. 



The following two lemmas are simple generalisations of Lemmas 2.1 and 2.3 



Lemma 3.2. There is a bijective correspondence between pyramids of size m > 1 and sequences 
(pi,P2, ■ ■ ■ ,Pr) of pyramids such that pi is a right Si-pyramid if i is odd and a left Si-pyramid if 
i is even, \p\ \ + • • • + \p r \ = m , s\ = 0, and 

1 < s i+i — Si < a — 1 if i is odd and 1 < Sj — Sj+i < a — 1 if i is even. 

Proof. If the pyramid p is not a right 0-pyramid there is a lowest piece in p above some interval 
that overlaps the negative real axis, that is an interval ]s2 — a, S2[, where < S2 < a. This piece 
is the bottom piece of a unique proper sub-pyramid p' and we can write p = p% p', where 
Pi is a right 0-pyramid. If p' is not a left S2-py r amid it contains a unique lowest piece above 
some interval ]s3, S3 + a[, where S2 — a < S3 < 82- This piece is the bottom piece of a proper 
sub-pyramid p" of p' , and we have p' = P2 Op", where P2 is a left S2-py ram id. The claim follows 
by repeating the argument a sufficient number of times. □ 

Lemma 3.3. There is a bijective correspondence between right 0-pyramids of pieces of length a 
and of size m and positive (am, m)- strings. 



Proof. The claim follows by a straight-forward generalisation of the proof of Lemma 2.3 the 



details of which are left to the reader. □ 

The following related correspondence between right pyramids and trees will not be needed 
in the proof of Theorem but may be of some independent interest. Recall that an n-ary 
tree, where n is a fixed positive integer, is a planar rooted tree all of whose vertices have order 
1 or n + 1 and whose root has order 1. The vertices of order n + 1 are called nodes. 
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Proposition 3.4. There is a bijective correspondence between right ^-pyramids of pieces of 
length a and of size m and a-ary trees with m nodes. 

Proof. First, note that there is an obvious bijective correspondence between (am, m)-strings 
and walks on the integers starting and ending at and consisting of m right-steps, each of 



length a — 1, and (a — l)m left-steps, each of length 1. In fact, t s as given in Definition 3.1 
defines the s'th site visited by the walk corresponding to a given string. Alternatively, the 
corresponding walk can be viewed as a path on the square lattice Z 2 from (0, 0) to (am, 0) with 
steps (l,o— 1) or (1,-1), called up-steps and down-steps, respectively. Positive strings then 
correspond to paths with vertices on or above the first axis only, and they are called generalised 
Dyck (a - 1) -paths [7], [8]. 

Consider a generalised Dyck (a — l)-path u and let uj' be the path obtained by removing 
the first step, which is necessarily an up-step. Thus uj' starts at height a — 1 and ends at height 
0. Let now uj\ be the part of uJ extending from the the first vertex in uj' at height to the 
final vertex (0, 0). Then uj\ is a generalised (a — 1)-Dyck path (possibly trivial), and uj' equals a 
path uj" starting at level a — 1, ending at level 1 and nowhere dropping below level 1, followed 
by first a down-step and then by uj\. Next, let UJ2 be the part of uj" extending from the first 
vertex in uj" at height 1 to the final vertex (x±,l). Then uj2 is a translated generalised Dyck 
(a — l)-path and the construction may be repeated a times to yield a decomposition of uj into 
a sequence uj a ,uj a -i, ■ ■ ■ ,uj\ of generalised (a — 1)-Dyck paths (suitably translated and possibly 
trivial) connected by single down-steps and preceded by an up-step. As a consequence, the 
number A m of generalised (o — 1)-Dyck paths with m up-steps satisfies the recursion relation 

A m = ^ A mi • . . . • A ma 

m-\ H \-m a —m — l 

m 1 ,...,m a >0 

Rephrased in terms of the generating function 

oo 

A(t)=Y,A m t m (3.1) 

m=l 

this relation takes the form 

A(t) = t(l + A(t)) a . (3.2) 

This identity has been noted previously in [7] (for a = 3). It is well known, and easy to 
establish, that the generating function for the number of a-ary trees as a function of the number 



of nodes likewise satisfies (3.2). We conclude that the coefficients are equal and hence, in view 



of Lemma 3.3, the claimed bijection is established. □ 



Corollary 3.5. The number A m of right pyramids of size m is given by 

1 / am\ (am)\ 

A m = - = — — - — — ^ m > 1. 

(a — l)m + l\ m J m!((a — l)m + 1)! 

Proof. It is known that the stated expression for A m equals the number of a-ary trees with m 



nodes, see e.g. [22]. Hence the claim follows from Proposition 3.4 □ 



Remark 3.6. It is, in fact, quite straight-forward to show directly, by a slight modification of 
the argument given in the proof of Proposition 1 in [6], that the generating function for the 



number of right 0- pyramids as a function of size satisfies the identity (3.2). The argument 
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involving generalized Dyck paths given above exhibits at the same time a proper generalisation 
of the well-known, and much exploited, correspondence between binary trees and standard Dyck 
paths, see e.g. [18]. 



In order to continue our efforts to establish a decomposition result analogous to Lemma 2.4 
some additional notation will be needed. We shall find it convenient to use the language of 
walks instead of strings in the following. Hence all walks subsequently will be assumed to have 
right-steps of length a — 1 and left-steps of length 1. A generic walk starting at i £ Z and 
ending at j £ Z will be denoted by Sij(m), S^Am!) etc. For i = the walk corresponding to an 



(n, m)-string is obtained by letting t s given as in Definition 3.1 be its s'th site. In Figures [2] and 



[3] and walks are illustrated by paths in Z 2 , replacing each right- step by an up-step (1, a — 1) and 
each left-step by a down-step (1, 1). In the case illustrated, a = 6. 

Given two walks S'^m') and S'- k {m"), the walk obtained by by first traversing Sy(m') and 
then S k Am") will be called the walk obtained by composing 5^(m') and S'- k {m") and will be 
denoted by S' i; j(m')S'! k (m"). Thus, composition of walks corresponds to juxtaposition of the 
corresponding strings. 

Evidently, an {am, m)-string is positive if and only if the corresponding walk takes place 
on the non-negative integers. Generally, we shall call a walk Sij{m) positive if j > i and the 
walk Sa{m) obtained by adding j — i left-steps at the end is a translate (by i) of a walk on the 
non-negative integers. Positive walks will be denoted by Pij{m),P' i -{m') etc. 

Given a walk Sij(m), its inverse walk Sjh^^Jm) is defined as the walk obtained by reflecting 
Sij{m) in the point j and reversing its direction of traversal. If i = j = this corresponds to 
reversing the order of the corresponding string. A walk is called negative if its inverse is a positive 
walk. Generic negative walks will be denoted by Nij{m),N' i -{m!) etc. Note that positive walks 
begin with a right-step whereas negative walks end with a right-step. 

For < j < i < a — 2 we denote by Tjj the straight walk from i to j consisting of i — j 
left-steps and, for 0<j<k<a— 2, we define 



Sij(m) = S' ik (m)T kj & S' ik (m) = Sij(m)U jk , 



that is the last equation means that Sij(m) contains at least k — j consecutive left-steps at 
the end, and S' ik (m) is obtained from Sij{m) by deleting its last k — j steps. In particular, we 
note that any walk Pa(m) necessarily ends with at least a — 1 left-steps such that Pa{m)Uik is 
well-defined, and ends with a left-step, for 0<i<k<a — 2. 

Definition 3.7. Consider a multiple composition of walks of the types Pa(m), Na(m), < i < 
a — 2, and Tij, 0<j<i<a — 2, and with possible insertions of terms Uik, 0<i<k<a — 2. 
By dropping the endpoint indices k and the step numbers m in the composition we obtain a 
word in the alphabet P, N, T, U. The composition is called admissible if only neighbouring pairs 
of letters of the form 

PN, NP, PT, TP, NT, TN, PU, UN (3.3) 
occur in the corresponding word. 

We are now in a position to formulate and prove the desired decomposition result for walks. 

Lemma 3.8. Any walk Soj(m), where < j < a — 2, can be written in a unique way as an 
admissible composition. 
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Proof. Let the walk Soj(m) be given and let us consider the two possible options, depending on 
the direction of its last step, separately. 



i) If Soj(m) ends with a right-step, any composition of the claimed type must be of the 
form 

S 0j (m) = S' Q j(m')Njj(m - m!) , 

where S' j(m') is either empty or is a composition that ends with a left-step, since walks corre- 
sponding to pairs in the list (3.3) not ending with N must end with a left-step. On the other 
hand, there evidently exists a unique S' j(m') ending with a left step such that the decomposition 
above holds. 



ii) If Soj(m) ends with a left-step, any composition of the claimed type must end with a 
Pjj(m'), a Tij or a Uij. 

If it ends with a Pjj{m!) and is not positive, the composition must have one of the following 
three forms: 

S Qj (m) = S'^m^Nj^m^Tj^Pj^^^j^Pj^m^) . . . l\ h l I ),{ ,),') , r > 1 , 

S Qj (m) = ,% r ( ///„ i A,. h ( m . i / > j:J . ( ///-j i'/), ; , /', ,;, { in-, i 7' , ; , / > jaja ( ///.-( i . . . T jrj P j:j (m') , r > 1 , 

S j(m) = S' oj (m )N jj (mi)Pjj(m') , 

where S' j(mo) is either empty or a composition of the claimed form, because P can only be 
preceded by T or N and T can only be preceded by P or N . Setting ji = j in the last case it is 
seen that in all three cases the last step in Nj 1 j 1 {m\) is the last (right) step, call it a, in Soj(m) 
whose initial point has negative value and whose final point is non-negative, and hence belongs 
to {0, 1, . . . , a — 2}. In case Soj(m) is positive we must have j = and Soj(m) = Pjj(m). 

That the P-and T-terms occurring in these compositions are uniquely determined can be 
seen as follows. Consider the step a defined above with endpoint j\. If the subsequent step is a 
left-step there exists a j2 < ji such that a is followed by Tj^ 2 and then by a right-step. This 
right-step is the initial step of a unique positive walk Pj 2 j2{ m 2) that is followed by a left-step, 
unless it equals Pjj(m'). If not, the argument can then be repeated. If the first step after a 
is a right-step, it is the initial step of a unique positive walk Pj 2 j2( m 2) that is followed by a 
left-step, unless it equals Pjj(m'). Now continue as previously until all P- and T-terms have 
been determined. 

In case the composition ends with a Tij or a Uij, essentially the same argument can be 
applied to establish uniqueness of the factors subsequent to the step a defined above. If a does 
not exist, i.e. if Soj(m) is positive, the unique composition of the claimed type must in this case 
be S j(m) = P o(m)U j. 

To establish existence of the composition for the part of Soj(m) subsequent to a one can 
proceed along the same lines just explained concerning uniqueness. Indeed, if a is followed by a 
left-step there must exist a non- negative ji < j\ such that a is followed by Tj 1 j 2 , which is then 
followed by a right-step. This right-step is the first step of a Pj 2 hk m 'i)- Choosing m<i maximal, 
it follows that Pj 2 n^ m 'i) i s either followed by a left-step, in which case the construction can be 
repeated, or the end of Pj 2 j 2 ( m 2) coincides with that of Soj(m), in which case the construction 
is finished, or the end of Pj 2 j 2 ( m 2) exceeds that of Soj(m) by j — 32 > left-steps, in which case 
these are annihilated by inserting Uj 2 j at the end. 

The case where a is followed by a right-step is treated in the same way. 

Together, i) and ii) prove the assertion of the lemma by induction. □ 
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Figure 2: Case i) of Lemma 3.8 String ending in N 

-v 





Figure 3: Case ii) of Lemma 3.8 String ending in PU 



Corollary 3.9. There is a bijective correspondence between walks of length am starting at 
with a right-step and ending at 0, and admissible compositions with initial term Poji( m i)> < 
ji < a — 2, mi > 1, and final term Poo(m r ), Noo(rn r ), m r > 1, or Tj r o, 1 < j r < a — 2, where 
mi + • • • + m r = m. 



4 Proof of Theorem 1.1 



The main purpose of this section is to give a proof of Theorem [Li] by exploiting the decomposition 
results of the preceding section. 

For m > 1, we let A m denote the number of positive walks Pu(m), which obviously is 
independent of i £ Z and also equals the number of negative walks Na{m). Moreover, A m 
also equals the number of right, respectively left, s-pyramids of size m as a consequence of 
Lemma 13.31 



By Lemma 3.2 it follows that the number B m of all 0-pyramids of size m can be written as 



(a-l) r - 1 A mi ...A r , 



(4.1) 



r>l WiH \-m r =m 



since the number of possible choices of the sequence (si, . . . , s r ) in Lemma 3.2 is (a — 1) 



\i — l 



On the other hand, it follows from Corollary |3.9| that the number of walks of length am 
starting at with a right-step and ending at can be written in the form 



r>l miH \-m r =m 



&r ^mi • • • -<4m r j 
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where r denotes the total number of P- and iV-terms, with sizes mi, . . . ,m r > 1, in a compo- 
sition and the factor a r counts the number of admissible compositions subject to the boundary 
conditions specified in Corollary 3.9 for fixed r and mi, . . . , m r . As indicated, this number only 
depends on r. Of course, the total number of walks of length am starting at with a right-step 
and ending at equals ( a ~-i) • Thus the proof of Theorem 
the following lemma. 

Lemma 4.1. For every a > 3 we have 



1.1 



will be completed by proving 



(o-l) 



r-l 



r > 1 . 



Proof. We use standard matrix techniques. Set b := a — 1 > 2 in the following and define the 
26 x 26-matrices £ , T and IA by 



£ = E < 



1 

1 



r 



i i 
i i 



u = u 



1 




where E is the b x b unit matrix, T is the lower triangular n x n-matrix with l's below the 
diagonal and O's elsewhere, and U is the transpose of T. We label the rows and columns of the 
first factor in the tensor products by i,j £ {0, 1, . . . , a — 2} and those of the second factor by 
R, S G {P, N}. Moreover, we shall use the ordering OP, ON, IP, IN,..., (a- 2)P, (a - 2)N of 
double indices, thus implying the standard identification of the tensor product of a b x 6-matrix 
and a 2 x 2-matrix with a 2b x 26-matrix. 

By construction, Smjs = 1 if and only if Rii{m)Sjj{m!) can occur in an admissible com- 
position, i.e. if i = j and R ^ S according to Definition 3.7 Similarly, TiRjs = 1 if and only 
if Rn(m)TijSjj(m') is allowed, and Uwjs = 1 if and only if Ra(m)UijSjj(m') is allowed. Viewing 
an admissible composition as a chain of links of one of the forms Ra{m)Sjj{m'), Rn(m)TijSjj(m') 
or Ra{m)UijSjj{m'), a chain of r — 1 links contains a total of r P- and iV-terms. Expanding 
the power (£ + T + U) r ~ 1 we hence get 

a r = {(£ + T + wr 1 } opop + {(£: + T + ^ 

3,S 

where the three terms correspond to the three possible types of final terms in the compositions 



specified in Corollary 3.9 Since Tjo = lifl<j<a — 2 and Too = 0, the result can be 
rewritten as 

/1\ 

a r = (1 0...0) {£ + T + U) r ~ l 

w 

Before using this to evaluate a r in general it is instructive first to consider the case a = 3 
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explicitly. In this case 
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such that 
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The characteristic polynomial of the 4 x 4-matrix entering this expression is found to be 

p 2 (A) = A(A-2)(A + 1) 2 

and hence its eigenvalues are 2, and —1. The eigenvalue multiplicities are seen to be 1 such 
that the matrix is not diagonalisable. 
In order to determine a r we define 



(a r \ 
a' r 

b r 



/o 1 o i\ 

10 

110 1 

\i i i oy 



r-l 



/A 

1 

1 

w 



which fulfills the recursion relations 

b r +l — Or+l 
K+l - br+l : 



a r 
b r 



b' 



i r+ i 



(4.2) 
(4.3) 



for all r > 1. Since b\ = b\ = 1 we get from the first equation in (4.3) that b r = b' r for r > 1. 
The first equation in (4.2) and the second in (4.3) imply 

b r = a r + a' r , r>2. (4.4) 

2b r , r > 2 , and 



Combining this with the last equations in (4.2) and (4.3) then gives 6 r +i 
hence 



h -2 



r-1 



3-2 



r-1 



Inserting this into the first equation in (4.2) yields 

O-r+l + Or = 3 • 2 r ^ 



r > 2. 



r > 1 . 
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The solution to this equation with initial condition a± = 1 is 



a r = 2 



r-1 



r > 1, 



as desired. 

We now return to the general case a > 3. Setting „4 = £ + T + W one can write A = X + y 
where X is the lower triangular 26 x 26-matrix whose matrix elements below the diagonal equal 
1 and are otherwise 0, and B is the upper triangular matrix with matrix elements equal to 1 in 
slots with even row and column indices above the diagonal and elsewhere, that is 

f0 1 1 1...1 1\ 

10 0...0 

110 10 1...1 1 

A= 1 1 1 0...0 

11111 1...1 1 
\l 1 1 1 1 1...1 1 0/ 

It is now easy to show by induction w.r.t. b that the characteristic polynomial of A is given 
by 

Ph (\) = \ b - 1 (\-b)(\ + l) b . 

Hence the eigenvalues of A are 0, b and —1, and it is readily seen that they all have eigenvalue 
multiplicity equal to 1. We denote by e the eigenvector with eigenvalue b normalised such that 
its first coordinate is b. One finds 

e = (b, 1, 6C, 2C, &C 2 , 3C 2 , ... , b( b - 2 , (b-l)Q h -\ b( b -\ bC'- 1 ), 

C = i + b-\ 



where 

Next, define 



fi :=(0, ... , 0, -i, 1, 1, ... ,1), l<i<6-l, (4.5) 

where — i is the 2{b — i)'th coordinate (such that the number of l's is 2i). Then f\ belongs to 
the kernel of A and one finds by direct computation that 

Afi = fi + f 2 + --- + fi-i, 2 < i < b - 1 . (4.6) 

It follows that the kernel of A h-1 is spanned by the vectors /i, ... , fb-i- By a straight- forward 
calculation one finds 



From this we conclude that 



(1, 1, 1, ... l)=b- 1 ^e-J2( b - 1 -' l f)j . 



b- 1 (100...0)X- 1 (e-^C &_1 "7^ 



= r 1 (l00...0)J _1 e 
= b- l b r - 1 (1 0...0) e 
= b r - 1 for r > 1, 
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where, in the second step, we have used (4.5) and (4.6) and the fact that fi has vanishing first 
coordinate for all i 

□ 



Recalling the definition of b the lemma is proven. 



Remark 4.2. As previously mentioned, the method used in this section to determine the number 
B m of pyramids of size m as given by Theorem |1.1| did not require knowing the number A m of 
right-pyramids of size m, as given by Corollary |3.5| For a = 2 the use of generating function 
techniques as in [6] proceeds by first determining the generating function A{t) for the A m and 
then using a simple algebraic relation between this function and the generating function 

oo 

B{t) = Y,B m t m 



for the B Tl 
form 



For general a > 2, this relation is a special case of eq. (|4.9|) below and takes the 
B{t) 



1 



(a-l)A(t) 



(4.7) 



For a = 2 the quadratic relation (3.2) satisfied by A(t) has a simple solution which easily yields 
the B m in closed form when inserted into (4.7). For a > 2 such a procedure does not seem 
feasible. 



We conclude this section by determining the asymptotic behaviour of the average width of 
pyramids of large size. For this purpose we first note that the asymptotic behaviour of 



B r . 



am — 1 
m — 1 



is readily obtained from Stirling's formula and is given by 

1 / a a \ m 

B m ~ — . r — - , m — ► oo . (4.8) 

y/27ra(a-l)m \(a - l)"" 1 J V ; 

Next, let B m ^ n denote the number of pyramids of size m and left width n, where the left width 
of a pyramid p with bottom piece covering the interval ]0,a[ equals n if the leftmost interval 
covered by a piece in p is ] — n, a — n[, and let 

B(t,v)= B n,mt m V n 
n>0,m>l 

be the corresponding generating function. Now recall the decomposition, in the proof of 
Lemma 3.2, of a pyramid p into a right 0-pyramid p\ and an arbitrary pyramid p' with bottom 
piece covering the interval ]s2 — a, S2[, where 1 < S2 < a — 1, or p' may be empty. This is seen 
to imply the relation 

B(t, v) = A(t)(l + {v + v 2 + --- + v a - l )B(t, v)) 
between A(t) and B(t,v), that is 

B i t ^) = - ? TTTTT- (4-9) 

v ; l-(v + v 2 + --- + v a - 1 )A(t) y ' 
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In particular, the relation (4.7) is obtained for v = 1. 

Moreover, differentiating eq. (4.9) with respect to v and setting u = lwe get 

C(t) = -a{a - l)B(t) 2 . (4.10) 
where C{t) is the generating function with coefficients 

Cm — ^ B m ,n • 
n>0 

The asymptotic behaviour of C m can now be obtained by standard singularity analysis. Indeed, 
from the polynomial equation (3.2) satisfied by A(t) we conclude that the singularity of A(t) 
closest to the origin is at 

_ (q - I)"-* 



and that 



where 



A(to) - A(t) ~ co 1 



i 

t \ 2 



to 



t — > to , 



co 



(a - l)- 2 v/2a(a- 1) and A(to) = (a - l)" 1 



Using (4.7) and (4.10) one obtains 

C{t) ■ 



1 / * 
4 V *0 



t — > to 



and hence [12] 



Cl7 



4*0 



1 



m — > cxd . 



(4.11) 



4 V(«- l) a_1 , 

For the average left width with respect to the uniform distribution of pyramids of size m we 
conclude from (4.8) and (4.11) that 



Br, 



-\/2ira(a — l)m , m 



oo . 



Hence we have proven the following result for the average width which equals twice the average 
left width plus a. 

Proposition 4.3. The average width of pyramids of size m is asymptotic to 



7T 



a(a — l)m 



for m — > oo. 
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5 Concluding remarks 



As noted in the introduction the pyramids under consideration in this paper may be considered 
as special planar LEGO structures built from 1 x a LEGO pieces. More specifically, consider 
the number of contiguous LEGO structures made out of m 1 x a pieces which are "flat" 
in the sense that all pieces are contained in the same vertical plane, and such that there is a 
unique piece in the lowest level of the structure. Pictorially speaking, the difference between 
that case and the one studied here is that pieces are allowed to hang underneath other pieces 
from the second level of the structure and upwards. These numbers turn out to be hard to 
compute; in [2] some of them have been calculated for 1 < a < 8 and are reproduced in The 
Online Encyclopaedia of Integer Sequences [21]. 

In [11] we have shown, using standard concatenation arguments, that the exponential growth 
rate 

lnL* 

g a = hm 

m— >oo Tfl 

is well defined and finite (and that this, in fact, holds also for more general non-planar classes 
of LEGO structures). Obviously, the asymptotic relation ( |4.8[ ) gives 

h a > ( ^ pT ~ e-Ha-l) (5.1) 

for large a, where we have set h a = e 9a . This lower bound appears to be rather tight; it is, indeed, 
at present our best lower bound for general a. It may be improvable by known techniques for 
fixed and relatively small values of a. For instance, for a = 2 one can appeal to the enumeration 
of multi-pyramids in [6] to get h<i > 9/2. By adapting the method of [11] to this setting we 
can improve this lower estimate further to /12 > 4.607 by computing the number c m of "fat" 
structures up to level 16 and proving c m +2 > 5c m for all m. 

Upper bounds on h a can be produced by adapting the method of [17] to this setting. Per- 
forming an analysis of depth 1 we can prove that, for all a, the largest root of 

(l a 9 _ y + f a 7 _ 3a 6 + 2a 5 _ 3 a 4 + l a 3 )2 ,5 + 
(_ 3a 8 + 17 a 7 _ 105 a 6 + 159 a 5 _ 73fl 4 + 165 a 3 _ 2^2 + 2a)x A + 

(_47 a 7 + 27a 6 - ifa 5 + fa 4 - ±f a 3 + |a 2 + \a)x*+ 

(a 6 - 4a 5 + 6a 4 - 4a 3 + a 2 )x 2 

is a majorant of h a , which in turn shows that 

h a < 6.356 a -4.375 

for large a. No closed form upper bound is available for analyses of depth 2 and 3, but ma- 
jorants are readily computable for a up to 8 as indicated Figure |5j Again, these appear to be 
approximately affine for large a. 

Using the Monte Carlo methods described in [1] and [2] we have produced estimates of h a 
for a up to 8. Strikingly, our estimates in each case have the form k a ( a _i)a-i ? with k a between 
1.238 and 1.264. This makes it tempting to speculate that 

5a a 

a ~ 4(0-1)--! ' 

In particular, it does not seem unlikely that hi = 5. Our best current estimate, achieved by a 
least square fitting of a function of the form AH n n c with Monte Carlo estimates for Lf e , . . . , L| 
yields H = 5.0012. 
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23456789 



Figure 4: Black: our conjectural formula for h a and the estimated values. Magenta: Upper 
bounds based on the Klarner-Rivest method (levels 1, 2 and 3). Blue: Lower bound from 
pyramids. Green: Lower bounds from counting fat buildings (levels 6 and 8). 
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